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ABSTRACT: This paper discusses research conducted by the U.S. Army Research Laboratory (ARL) - Vehicle 
Technology Directorate (VTD) on advanced suspension control. ARL-VTD has conducted research on advanced 
suspension systems that will reduce the chassis vibration of ground vehicles while maintaining tire contact with 
the road surface. The purpose of this research is to reduce vibration-induced fatigue to the Warfighter as well 
as to improve the target aiming precision in-theater. The objective of this paper was to explore the performance 
effectiveness of various formulations of the Generalized Predictive Control ( GPC) algorithm in a simulation 
environment. Each version of the control algorithm was applied to an identical model subjected to the same 
ground disturbance input and compared to a baseline passive suspension system. The control algorithms 
considered include a GPC with Implicit Disturbances, GPC with Explicit Disturbances, and GPC with Preview 
Control. A two-axle tactical vehicle with independent front and rear suspensions was modeled in the TruckSim 
full-vehicle dynamics simulator. The control algorithms were compared based on their effectiveness in 
controlling peak acceleration and overall average acceleration over a range of vehicle speeds. The algorithms 
demonstrated significant reductions in the chassis acceleration and pitch of the full -vehicle model. 
Keywords: Vehicle vibration control, vehicle dynamics simulation, control algorithm, advanced suspension, 
generalized predictive control. 



I. INTRODUCTION 

Passive suspension systems for Army wheeled vehicles are optimized for off -road passenger comfort 
and maneuverability, while maintaining acceptable on-road handling. These two goals often lead to conflicting 
design requirements. With a view toward available suspension systems that would improve both performance 
metrics, the Vehicle Technology Directorate (VTD) of the U.S. Army Research Laboratory (ARL) has 
underway a research program aimed at developing active suspension systems for implementation on the next 
generation of Army vehicles. The overall objective is to reduce vibration -induced fatigue as well as improve 
target aiming precision in-theater, while maintaining on-road handling. There are basically two main groups of 
active suspension systems: semi-active and fully active. Semi-active suspension systems vary spring and damper 
properties to store or remove energy from the system according to some applied control algorithm. Usually, this 
is achieved through variable damping components such as variable orifice valves or magneto -rheological fluids. 
Fully active suspension systems use force actuators to add and remove energy from the system. This can be 
achieved using hydraulic, pneumatic, or electromechanical actuators. Current research at VTD is addressing 
both types of systems. However, this paper focuses on fully-active controls and the application of a predictive 
control technique known as Generalized Predictive Control (GPC). GPC is a linear, time -invariant, multi- 
input/multi-output method that uses an Auto Regressive exogenous (ARX) model to represent the input-output 
behavior of the system and to design the controller [3, 4, 8, 7, 9]. Three variants of the techniques have been 
evaluated. In the first variant, no explicit account is taken of the disturbances acting on the system during system 
identification and controller design [4, 8]. This version is called GPC with Implicit Disturbance in this paper. 
The second and third variants take explicit account of measurable disturbances. The first of these two (called 
GPC with Explicit Disturbance herein) uses only the past values of the measured disturbances while the latter 
(called GPC with Preview in this paper) additionally uses measurements (or predictions) of future disturbances 
[3, 7, 9]. The purpose of this paper is to present the results of numerical simulations of the application of the 
aforementioned GPC variants to a multi-body dynamic model of an Army tactical vehicle with nonlinear 
suspension systems. The paper begins with a summary of the key equations describing the underlying method 
including system identification, computation of the control law, and calculation of the commands sent to the 
actuators. Then, a brief discussion of the simulation environment and the vehicle model are presented. This is 
followed by a description of the three road profiles used as disturbances in the simulations. Representative 
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results from the numerical simulations are then presented. The paper finishes with some concluding remarks on 
the results of the numerical study and comments on the nature of some planned future work. 



II. CONTROL METHODOLOGY 

The steps in developing a GPC-based controller are system identification and derivation of the control 
law. System identification is based on the derivation of a multi-step output prediction equation using the 
Observer Markov Parameters (OMP) that comprise the ARX model that is used to characterize the system. The 
control law then follows from the minimization of a quadratic performance equation involving the predicted 
response and the control command. A summary of the key steps of the process, adapted from the more extensive 
development present in reference [4], is given below. Since GPC is a multi -input/multi-output controller, the 
control command, fc f }, measured response, [»}, and external disturbance, [d t } can be vectors whose elements 
correspond to the values of a specific actuator or sensor at time t. 
GPC with Implicit Disturbance 

The relationship between the input and output time histories of a linear, time -invariant multi-input/multi-output 
system can be described by an ARX model of the form: 
y(t) = a iy (t - 1) + a 2 y(t - 2) + ... + cc p y(t - p) + 

y? D «(0 + ftufe - 1) + ... + p p u(t - p} 

This equation says that the output, y(t), at the current time step, t, may be estimated by using p sets of the 
previous input and output measurements, u(t-l), u(t-p), y(t-l), ...,y(t-p), and the current input, u(t). The 
integer p denotes the order of the ARX model. The coefficient matrices Di and Di are called Observer Markov 
Parameters and are the quantities computed by the identification process [2]. 

To start the identification process, the system is excited with band-limited white noise. These independent 
random excitations are applied concurrently to all of the r c control inputs. The r c control inputs and the m 
response outputs are recorded for L samples. The resulting input and output time histories, u and y, are then used 
to form the data matrices, y and V. 
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[r c + (r e + m)p] x L 

Throughout this paper, the dimensions of the matrices and vectors are noted below key equations. The order of 
the ARX model, p, and the number of samples, L, are specified by the user. Some guidelines for their selection 
are given in reference [2]. The matrix V, containing the OMP, Dj and Dj, follows from equation 2 by inversion 
and has the form: 

F=[£ D fit ff! p 2 cr 2 ... £ p 

mx [r e + O ff + m) p] 

The one-step-ahead output prediction equation given by equation 1 is the starting point for deriving the multi- 
step output prediction equation that is needed for deriving the GPC controller. This equation is obtained by 



ffp] 



(5) 



replacing t by t+j in equation 1, and letting j range over the set of values /=/, 2, ... h p -l, where h p is the 
prediction horizon (the number of time steps for which the future responses are predicted). The resulting 
equations can be assembled into a multi-step output prediction equation having the form: 
(yjj = + a(y p ) + p{«p) (6) 

The coefficient matrices □ , □ , and Dare formed from combinations of the OMP. The quantity >^ p is the vector 
containing the future predicted responses, whereas -u^ is the vector containing the (unknown) future control 
commands. The integer h c is the control horizon (the number of time steps over which future control is assumed 
to act). The quantities u p and y p are vectors containing the previous p sets of control commands, and measured 
responses, respectively. 
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where the coefficients are calculated using the following relations: 

In this formulation the size of the ARX model, /?, is set equal to the number of future predicted responses, h p , 
(prediction horizon) and future commands, h c , (control horizon). Generally, these three parameters can differ 
from one another [3, 8]. 

The optimal control law is obtained by minimizing an objective function. To accomplish this, one defines a 
quadratic objective function as follows, 

/= y^Qy fep +<^ (9) 

where 

Q is a symmetric positive semi-definite matrix that assigns weights to the predicted responses, 
R is a symmetric positive definite matrix that assigns weights on the future control commands, 
The superscript T indicates the vector or matrix transpose. 
Substituting the expanded equation for [y] yields, 

The minimum is obtained by taking the partial derivative of the function with respect to future commands, {u^} 
and setting it equal to zero. 

0 = ^-[xf« A J+ *{y P }+ P{«p}]' Q[r{u hr ) + a{y p } + pf^}] + ^h^JIh^ 

(lib) 



(10) 



(11a) 



0 = 2x T Q [t{«,J + a[y p } + p{«p}] + 2JI j : 

Solving for control commands, yields the following optimal future commands: 
{^^-(^QT+fO-^Qiaiy^+^Up)], " (12) 
Defining the following control law matrices, 

A c = -(-T^T-h Rj-VQfl (13a) 
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b c = -(t t Ot + r]-Vgj? (13b) 

yields control laws that calculate the next h c sets of commands, 

W = A £ ty p }+Bj« p } (14) 

where corresponds to the control commands h c time steps into the future. 

However, only the first set of control commands will be applied to the system, thus only the first r c rows of 
[uft } are needed. The remaining rows are discarded. 

r«i-i] 
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GPC with Explicit Disturbance 

Generalized Predictive Controller with Explicit Disturbance, as discussed in [3, 9], is similar to GPC with 
Implicit Disturbance, except that now the external disturbances are explicitly measured and included in the 
model. Time histories of controls commands, measurement responses, and, now, disturbance measurements are 
used to generate the ARX model for the system. The V matrix with disturbance information is shown below. 
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lTd + r c + f r d + r c + ™')p] X L 
The new matrix V , containing the OMP has the form: 



(17) 
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This OMP values are used to predict future responses based on past responses (including external disturbances) 
and control commands. 



b r t + fc p J l^t+jj (yt-pJ i«r-pJ 




(18a) 



(18b) 



As before, the next control commands are calculated by minimizing the objective function and taking the partial 
derivative with respect to control commands; the rows corresponding to future control commands are dropped. 
Using this extension, the measured disturbance information can be used. 
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GPC with Preview Control 

Preview Control or Look-Ahead Control measures the external disturbance before it affects the vehicle and 
incorporates this information into the control algorithms. Preview Control was first published in [1], and is 
further developed in [5, 6] for linear full-state feedback control. 

Generalized Predictive Controller with Preview Control is a direct extension of GPC with Explicit Disturbance. 
If it is possible to generate a model of the explicit disturbance in order to predict future disturbances [9]; this 
future information can be used in conjunction with the control command, external disturbance, and measured 
response time histories to generate the future control commands. This research, however, assumes that the future 
disturbance can be directly measured using vehicle based sensors. Again, time histories of controls commands, 
measurement responses, and disturbance measurements are used to generate the ARX model for the system. The 
matrix Y is identical to equation (17). These are used to predict future responses based on past responses, past 
and future disturbances, and past and future control commands. In this formulation, it is assumed that the current 
and future disturbance can be measured. The matrix V is the component of the ARX model that corresponds to 
the current and future disturbances. 
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As before, the future control commands are calculated by minimizing the objective function and taking the 
partial derivative with respect to control commands; only the rows corresponding to the next immediate control 
commands are retained. Using this extension, both the future and the previous measured disturbance information 
can be used. 

fri f+ o) (Yt-i) pvi] K-0 

[*r + o}=y c i Uofcj : +fcj : +6 C = (21) 
{dt + hj bt-pJ U*r- p J (dt-pj 

The difference in control command between GPC with Explicit Disturbance and GPC with Preview Control is 
solely due to the presence of the future disturbance term. The control matrices, D c , U c , and U c are the same for 
both control laws. 
Modeling and Simulation 

TruckSim is a widely utilized, commercially available software tool used to simulate and analyze the dynamic 
behavior of wheeled vehicles. It is employed in conjunction with Matlab/Simulink in order to simulate the 
effects of the GPC algorithms on a full- vehicle model. 




Figure 1 - TruckSim Modified HMMWV model on Ferryman 3. 



The Truck Sim model takes the full vehicle and terrain into consideration, including: driver controls and 
steering, road geometry, suspensions, tires, and powertrain. Component parameters are input into TruckSim via 
constants, linear coefficients, non-linear tables and algebraic formulas. 

The driver model controls the vehicle steering and throttle in order to follow a specified path and speed profile. 
Truck Sim removes artificial vehicle constraints; therefore, vehicle speed and direction can deviate from the 
specified profile. 

Full three-dimensional (3-D) road geometry is specified at 0.1 meter increments. The steering system captures 
the steer of each wheel due to the steering system geometry and compliances. The suspension system models 
capture the full non-linear kinematical behaviors and compliances of the front and rear suspensions. These 
models include data for springs, dampers, and jounce/rebound stops. The tire model is an internal table -based, 
single point contact model. Non-linear tables using actual measured tire data represent vertical force, lateral 
force, longitudinal force, aligning moment, and overturning moment as functions of deflection, slip, load, and 
camber. The powertrain model includes engine torque, torque convertor characteristics, and transmission gear 
ratios and efficiencies. A representative TruckSim tactical vehicle model is shown in Figure 1. 
TruckSim provides a generalized two -axle tactical vehicle model (HMMWV) with an independent front and 
rear suspension system. Figure 2 contains a graphical layout of the suspension as modeled in SuspensionSim. 




Figure 2 - Representative HMMWV Suspension 
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For this research, we modified some of the basic HMMWV vehicle parameters to study the performance of GPC 
algorithms on various terrain profiles. The vehicle mass was doubled. The suspension springs were stiffened, 
and the front damping coefficient was reduced. This is shown in Table 1 . 



Table 1 - Vehicle parameters for modified HMMWV model 





Basic 
HMMWV 


Modified 
HMMWV 


Sprung Mass 


2210 kg 


5000 kg 


Front spring 


150 N/mm 


250 N/mm 


Front damper 


30 kNs/m 


15 kNs/m 


Rear spring 


150 N/mm 


400 N/mm 


Rear damper 


30 kNs/m 


30 kNs/m 



The equations of motion in the TruckSim math models are ordinary differential equations (ODE) 
derived from the first principles for full nonlinear 3-D motions of multiple connected rigid bodies. The sprung 
mass of the vehicle is a rigid body with six degrees of freedom (DOF). Each wheel has one spin DOF. Each 
independent suspension axle has two multi-body DOF for the vertical movements of the wheel centers. Other 
suspension motions such as camber, lateral position, etc., are constrained as nonlinear functions of the 
independent variables. Each suspension has six compliance DOF. Each tire has two dynamic DOF, one for 
lagged lateral slip and one for lagged longitudinal slip. 

Because the GPC algorithms are implemented in Matlab and Simulink, TruckSim uses co-simulation to update 
the math model and the algorithm calculations simultaneously. TruckSim provides an S -function block to 
communicate with Simulink. The control algorithms run in Simulink using the discrete integration method with 
a fixed time-step of 0.001 seconds (1000 Hertz). TruckSim uses the Adams -Moulton second-order integration 
(AM-2) method; the time-step for this integration method is fixed at half the Simulink simulation rate. The 
model is integrated before and after each communication step with Simulink. 

Within Simulink, the GPC algorithms sample the vehicle sensors once every 0.01 seconds (100 Hz). The 
actuator commands to TruckSim are, also, updated every 0.01 seconds. For this work, the control algorithms 
optimize commands for force actuators at each of the four suspension corners to minimize the sprung mass 
motion. TruckSim uses the External Springs option to introduce outside forces from the actuators. Additionally, 
the absolute road heights at the tire contact patch and at the preview sensor location are sent to the algorithms. 

Road Profiles : 

Three road profiles were used to evaluate the algorithms: Sine, Curb, and Perryman 3. In each case, the 
road profiles were provided as an array of ground distances and vertical heights. 
Sine 

This road profile was defined by a sine wave with a half-amplitude of 0.10 meters and a wavelength of 10 
meters. The intent was to provide a single -frequency disturbance that relates to vehicle speed, in this case, 1 
meter/second (3.6 kph) of speed corresponds to 0.1 Hz excitation. There is a 30 meter flat stretch to allow the 
model to reach steady state after startup. Figure 3 shows a portion of this type of road profile. 



Road Distance - m 



Figure 3 - Sine road profile. 

Curb 

This road profile is defined by a 0.10 meter step rise followed by 10 meters of level surface, and then a 0.10 
meter step drop, as shown in Figure 4. The intent was to provide a step input disturbance into the system. There 
is a 30 meter flat stretch to let the initial oscillations die out and the vehicle reach steady state. 
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E 01 



Road Distance - m 



Figure 4 - Curb road profile. 

Perryman3 

Perryman 3 is an off -road test track located at the Aberdeen Proving Ground. This road profile, shown in Figure 
5, is intended to evaluate the algorithms against a challenging, "real-world" course. 




Road Distance - m 



Figure 5 - Perryman 3 road profile. 

Simulation Results 

Simulations were run until the vehicle traveled 150 meters. The open loop model was run in parallel with 
closed-loop GPC with Implicit Disturbance, GPC with Explicit Disturbance, and GPC with Preview Control for 
performance comparison studies. All measurements assumed no noise or offset-bias. Initial conditions for the 
suspension system state variables were set to zero. Vehicle speed was held at a constant value for the entire run. 
For this work the system identification for each algorithm was pre-computed. Band-limited white noise was 
used to excite the four actuators. The vehicle was run over Perryman 3 terrain profile at 40 kph. This provided 
broadband random disturbance inputs out to 50 Hertz with large peaks from 0.5 to 3 Hertz. 
The following objective functions were used as metrics in several comparisons of algorithm performance: 

= lowpass(^y r y } (22a) 
Actuator Force = low pass(>fu T u~) (22b) 
Actuator Power = low pass^lp ower J power ] (22c) 

where "lowpass" denotes a 4th order digital low-pass Butterworth filter with a cutoff frequency of 1 Hz. 
Using the cosine double angle formula, it can be shown that the square of a signal will have a steady offset 
component and a high frequency component that is twice the frequency of the original signal. The low-pass 
filter removes the oscillatory component, leaving only the steady offset, which is related to the amplitude of the 
fast-fourier transform of the original signal. Because the cutoff frequency is higher than the lowest frequency of 
excitation in the road profiles, there is a small low frequency oscillation in the objective function, J. The low 
frequency oscillation is removed by averaging the objective function over the entire run-time. 
For each simulation, the averaged values of J were plotted for a series of vehicle speeds up to 80 kph. The 
closed-loop responses were normalized by the corresponding values of the open-loop responses at each speed. 

Response Improvement = 1 — (23) 

/iC-esp. ppen [npp 

This normalization allows a more meaningful comparison of the algorithms despite the variation in actuator 
force or power. 

Force Normalized _ RssvffnseJTKprovemffrit^Tritj^ 

Effec t iven ess Actuator Fore tui^ritk-m 

Power Normalized Resutmse iTtzprovemejit^^q^ 

Effectiveness 



(24a) 
(24a) 



Actuator Pemsr Rtg „ lAm 

Simulation results based on the Modified HMMWV parameters are presented in Figures 6 through 12. These 
figures were generated using the following algorithm parameters, as shown in Table 2. 
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Table 2 - Algorithm tu ning parameters used for modified HMMWV vehicle. 





GPC with 

Implicit 

Disturbance 


GPC with 

Explicit 

Disturbance 


GPC with 

Preview 

Control 


actuators, r c 


4 


disturbances, r d 


4 


sensors, m 


3 


Q 


50 (vertical acceleration), 10 (pitch), 1 
(roll) 


R 


range { le-7 to 5e-6 } 




ARX model size, p 


8 


Prediction 
Horizon, h p 


8 


Control Horizon, 

h c 


8 


Preview Horizon 






8 (0.08s) 


ID sample size, L 


400 points 




Disturbance for ID 


Perryman 3 at 40 kph 


Actuator ID signal 


white noise 


Actuator Limit 


1.5e5N 



Effects of Control Weighting 

The modified HMMWV model was run with GPC with Preview Control over the Sine road profile from 
0 to 80 kph. The objective function, actuator force, and actuator power are averaged over a 150 -meter 
simulation. The speed range corresponds to a sinusoidal excitation of 0 to 2.2 Hz. This is repeated several times 
for a range of control weights (R matrix values); the response weights (Q matrix values) were held constant. 
This effect of the Control Weighting matrix is shown in Figures 6(a) through 6(f). 

An order of magnitude decrease in the control weight matrix produces a significant improvement in response 
and a corresponding increase in the required force and power. The force normalized effectiveness shows that the 
curves for each weight coalesce in the same regions; the exception of R = 0.1e-6 occurs because the forces 
required are at the limit for the actuator. This suggests that as long as the commanded forces are within the 
capability of the actuator, the level of the response is linear with control loads. However, when comparing the 
actuator power, the power-normalized effectiveness decreases as the control weights decreases. This shows that 
the actuator power costs increase much faster than the responses decrease. Actuator power, then, becomes the 
driving factor for choosing control weights. 
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Vehicle Speed - kph 



R= 5.0e-6 
R= 1.0e-6 
R=0.5e-6 
R=0.1e-6 




Vehicle Speed - kph 



(e) if) 

Figure 6 - Effect of control weights on GPC with Preview Control on sine road profile. 

Comparison of Algorithms 

A short comparison of the three control laws operating on the three road courses (disturbance inputs) is 
presented in the discussion below. The control weight for GPC with Explicit Disturbance was R=1.0e-7 over the 
entire speed range. To simplify analysis, the control weights for the other two algorithms were chosen such that 
the magnitude of the measured responses matched GPC with Explicit Disturbance at each speed. Again, the 
objective function and the actuator force are averaged over the 150-meter run for sinusoidal disturbances from 0 
to 80 kph (corresponding to 0 to 2.2 Hz). Figures 7(a) through 7(d) show a comparison of the three algorithms 
versus vehicle speed. 

The open loop response reaches a peak at 50 kph; representing the vehicle's pitch mode resonance 
frequency. The closed-loop algorithms reduce the peak by 60%; effectively attenuating the natural frequency. 
The normalized data shows that, GPC with Explicit Disturbance and GPC with Preview Control have similar 
performance; requiring the same actuator force and power. However, GPC with Implicit requires significantly 
greater actuator force and power, almost twice the actuator power at 80 kph. The addition of road disturbance 
information to the GPC algorithm greatly reduces the actuator loads for the closed loop system. 
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(c) (d) 
Figure 7 - Effect of GPC algorithms on sine road profile. 



Figures 8(a) through 8(d) show the time-history of the measured responses and actuator forces for each 
algorithm as the vehicle maneuvers over the Sine road profile at 50 kph, which represents a single frequency 
disturbance of 1.4 Hz. 

All three versions of the GPC algorithm reduce the chassis vertical acceleration and pitch. Chassis roll increases, 
due to the very low response weight assigned to chassis roll. GPC with Implicit Disturbance performs worse on 
controlling chassis vertical acceleration, but slightly better at reducing chassis pitch. Since the chassis pitch is 
the dominant term in the objective function, this results in similar levels for the response objective function. It is 
clear that GPC with Implicit Disturbance requires much higher actuator loads than the other two algorithms. 
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(c) (d) 

Figure 8 - Time History of Control Algorithms for Light Vehicle parameters on Sine course at 50 kph. 



The same simulations were run using the Perryman 3 road profile; the algorithm performance is shown in 
Figures 9(a) through 9(d). The GPC algorithms reduce the overall response by more than 25% from the highest 
open loop response. Much like the performance on the Sine road profile, GPC with Explicit Disturbance and 
GPC with Preview Control have almost identical response and actuator objective functions. Once again, GPC 
with Implicit Disturbance requires significantly higher actuator forces and loads to match the response objective 
function. 
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Figure 9 - Effect of GPC algorithms for Modified HMMWV on Perryman 3 road profile 



The GPC algorithms reduce the peaks of the chassis pitch compared to open loop cases. This dominant sensor 
makes up for the slight increases in closed-loop roll and vertical acceleration. Without measured knowledge of 
the disturbance, GPC with Implicit Disturbance must wait until the sign of the responses change before applying 
a countering force. This results in the large, but brief forces as the algorithm catches up to the changing slope of 
the disturbance. The algorithms with disturbance information have a much quicker response, thus reducing the 
peaks. The time history of the vehicle on Perryman 3 is shown in Figures 10(a) through 10(d). 
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The performance of the algorithms for the Curb road profile is shown in Figures 11(a) through 11(d). Again, all 
three GPC algorithms reduce the response objective function, by as much as 14% at the highest speed. However, 
this time, there is a surprising difference between the actuator force and power of the algorithms. GPC with 
Explicit Disturbance requires the highest actuator loads, double the force and power as required for GPC with 
Implicit Disturbance. GPC with Preview Control requires the least power of the three variants; less than 200 
Watts above 50 kph as shown in Figure 1 1(d). 
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Figure 11 - Effect of GPC algorithms on curb road profile. 



Previously, GPC with Explicit Disturbance was nearly indistinguishable from GPC with Preview Control. 
Figures 12(a) through 12(d) provide a closer look at the actuator forces experienced as the vehicle encounters 
the step height increase of the road. When the GPC with Explicit Disturbance encounters a step increase in the 
disturbance height (basically, a very steep slope), this steep slope produces a high control command. The steep 
slope is in the disturbance history for 8-timesteps. On the other hand, GPC with Preview Control issues small 
commands (100 N) for the 8-timesteps prior to the start of the curb. This reduces the peak force required once 
the tire contacts the curb, resulting in the lowest peak force. Additionally, after the tire contacts the curb, the 
future disturbance contains 8 time-steps of flat data to further temper the control commands. 
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As before, GPC with Implicit Disturbance is unaware of the disturbance and can only respond after the 
measured responses are affected by the curb. Thus, it commands high peaks, but not in response to an infinite 
slope. 
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Figure 12 - Time History of GPC Algorithms for Light Vehicle parameters on Curb course at 72 kph. 



Concluding Remarks 

The results of numerical simulations directed at evaluating the effectiveness of a predictive control 
method known as Generalized Predictive Control (GPC) for active control of chassis suspension systems in 
Army wheeled vehicles have been summarized. GPC is a linear time -invariant, multi-input/multi-output 
predictive control method that uses an ARX model to describe the input-output relationship characterizing the 
system. The coefficient matrices of the ARX equation are determined using system identification techniques and 
then used to form the matrices comprising the GPC control law. Key equations in the formulation of the method 
were summarized. Three variants of the basic GPC algorithm were investigated. The first assumes that the 
disturbances are unknown and computes the control law taking implicit account of the disturbances acting on 
the system (GPC with Implicit Disturbance). The second assumes that the past values of the disturbances can be 
measured and uses that information in deriving the control law (GPC with External Disturbance). The third 
variant assumes that the disturbances ahead of the vehicle can be measured and additionally uses that data when 
deriving the control law (GPC with Preview Control). Closed-loop behavior of the system was simulated using a 
full- vehicle nonlinear model of a two -axle tactical vehicle moving over three types of road profile disturbances 
for each of the three variants of GPC and compared to the open-loop behavior of the system. These comparisons 
show significant reductions in chassis motions of the model. Through simulation, it has been shown that, for a 
given algorithm, the control weights reduce sensor objective function at the same rate that it increases actuator 
objective function. Thus the control weight can be directly used to scale the total force available to a controller. 
GPC with Implicit Disturbance consistently requires higher actuator forces and more power than GPC with 
Preview Control or GPC with Explicit Disturbance when given a tonal or broad band disturbance. GPC with 
Preview Control performs nearly identically to baseline GPC with Explicit Disturbance; however, the preview 
control algorithm does drastically reduce the initial pulse due to step inputs. Despite the presence of disturbance 
information, GPC with Explicit Disturbance performs poorly on step road profile, even worse than having no 
disturbance information at all. In summary, for the three road course profiles investigated, GPC with Preview 
Control performs better overall than the other GPC variants for reducing vehicle -chassis acceleration and pitch. 
The simulation results show that GPC with Preview Control is capable of reducing the vehicle-chassis 
acceleration by as much as 25% on average for the most challenging terrain profile course simulated. The 
differences in the performance of the various algorithms studied in this research effort are currently being 
examined further to understand the main causes for the observed performance differences. In future, these GPC 
algorithms will be tested using a two-degree-of freedom suspension test rig at the Army Vehicle Research 
Laboratory to verify and validate these algorithms. 
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ARL 

ARX 

d 
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GPC 
h p , h c 
Hz 
kph 
L 

LVDT 
m 

m, r c , r d 
N 

OMP 
P 

Q,R 
s 
t 
u 

u(t), y(t), d(t) 
u(t-j) ,y(t-j), d(t-j) 
u(t+j) ,y(t+j), d(t+j) 
V 

VTD 

y 

Y 

□i, □ 

□ C D □ QD □ QU □ Q 

□ □□□□□□□□□ 



U.S. Army Research Laboratory 
AutoRegressive exogenous input 
disturbance inputs 
control actuator force 
Generalized Predictive Control 
prediction and control horizons 
Hertz 

kilometers per hour 

number of samples used for system 

identification 

Linear Velocity Displacement 

Transducer 

meter 

number of response outputs, control 
inputs, disturbance inputs 
Newton 

Observer Markov Parameters 
order of the ARX model 
response penalty, control penalty 
seconds 

index of current time step 
control inputs 

vectors of current inputs, outputs, 
and disturbances at time step t 
vectors of past inputs, outputs, and 
disturbances at time step (t-j) 
vectors of future inputs, outputs, and 
disturbances at time step (t+j) 
system identification data matrix 
Vehicle Technology Directorate 
response outputs 

matrix of observer Markov 
parameters 

observer Markov parameters 
control law gain matrices 
coefficient matrices in multi-step 
output prediction equation 
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